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Abstract 

The accelerating expansion of the Universe at the present stage is a process that will 
change the spectrum of relic gravitational waves. Here we present a numerical calculation 
for the power spectrum of relic gravitational waves in the accelerating Universe. The 
results show that although the overall features of the power spectrum is similar to that 
in the non-accelerating models, the amplitude is smaller by an order of 1CU 1 . We also 
find that the spectrum is very sensitive to the index f3 of the inflationary expansion with 
the scale factor o(r) oc |t| 1+ ^. With increase of /3, the resulting spectrum is tilted to be 
flatter with more power on high frequencies, and the sensitivity of the second science run 
of the LIGO detectors puts a restriction on the parameter /3 < —1.8. The influence of 
reheating following the inflation has been examined. 

PACS numbers: 98.80.-k, 98.80.Es, 04.30.-w, 04.62+v, 

Key words: gravitational waves, accelerating universe, dark energy 

e-mail: yzh@ustc.edu.cn 



1 



In the past the stochastic background of relic gravitational waves has been extensively studied [1] 
[2] [3]. The details of the spectrum of relic gravitational waves depends not only on the early stage 
of inflationary expansion, but also on the expansion behavior of the subsequent stages, including 
the current accelerating expansion. The calculations of spectrum so far [4] [5] [6] [7] [8] [9] have 
been done for the transitions from the inflationary era to the radiation, or matter dominated state, 
both being a decelerating expansion. However, the recent observations [10] [11] [12] indicate that the 
Universe is currently under accelerating expansion, which may be driven by the cosmic dark energy 
(Q\ ~ 0.7) plus the dark matter (Cl m ~ 0.3) with Q = 1 [13]. The current accelerating expansion 
of Universe will have an impact on the relic gravitational waves. We have studied this effect and 
obtained an approximate spectrum of the relic gravitational waves using an analytic method [14]. 
In this paper we present a numerical calculation of the spectrum. The resulting spectrum from the 
numerical calculation agrees with that from the analytic one. Throughout the paper we use units 
with c = h = 1 and adopt notations similar to that of Grishchuk [6] for convenience for comparison. 

The history of overall expansion of the Universe can be modelled as the following sequence of 
successive epochs of power-law expansion. 

The initial stage (inflationary) 

o(r) = l | r -oo < t < n, (1) 

where r is the conformal time, 1 + f3 < 0, and t\ < 0. The special case of j3 = — 2 is the de Sitter 
expansion. 

The z-stage 

o(r) = a z (r - Tp ) 1+f> ' , n<T<r s . (2) 

The reheating process towards the end of inflation can be quite complicated and is model-dependent. 
This z-stage for various f3 s can represent a general reheating epoch . 
The radiation-dominated stage 

a(r) = a e (r - r e ), r s < r < r 2 . (3) 

The matter-dominated stage 

a(r) = a m (r - r m ) 2 , r 2 < r < t e , (4) 

where te is the time when the dark energy density p\ is equal to the matter energy density p m . 
Before the discovery of the accelerating expansion of the Universe, the current expansion was usually 
taken to be in this matter-dominated stage, which is a decelerating one. The value of the redshift 
ze at the time te is given by (1 + z E ) = o,{th)/ol{t e)i where th is the present time. Since p\ is 
constant and p m (r) oc a~ 3 (r), one has 1 = — = — -, — — If the current values ~ 0.7 
and $7 m ~ 0.3 are taken, then it follows that 1 + z E = (^f 1 ) 1 ^ 3 ~ 1-33. 
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The accelerating stage (up to the present) 

a(r) = l H (T a - ry 1 , t e <t<t H - (5) 

This stage describes the accelerating expansion of the Universe, which is the new feature in our 
model and will induce some modifications to the spectrum of the relic gravitational waves. It should 
be mentioned that Eq.(5) is an approximation to the current expansion behavior since the matter 
component also exists in the current Universe. 

Given a(r) for the various epochs, the derivative a' = da/dr and the ratio a' /a follow immediately. 
There are eleven constants in the above expressions of a(r), among which the /3 is imposed upon as the 
inflationary model parameter. By the continuity conditions of a(r) and a(r)' at the four given joining 
points n, t s , T2, and te, one can fix only eight constants. Besides, choosing the overall normalization 
of a at the present time th and taking the expansion rate as the observed Hubble constant will fix 
two other constants. Specifically, we put (r a — th) = 1 as the normalization of a(r) at the present 
time th, which fixes the constant r , and the constant Ih is fixed by the following calculation 

so Ih turns out to be just the Hubble radius at present. Thus after a value of the model parameter 
/3 is picked up, all the eleven constants in a(r) are fixed up. In the expanding Robertson- Walker 
spacetime the physical wave length is related to the comoving wave number by 

A S ^, (T) 
and the wave number kn corresponding to the present Hubble radius is 

k H = = 27r . (8) 

Incorporating the perturbations to the spatially flat Robertson- Walker space-time, the metric is 

ds 2 = a 2 {r)[dT 2 - (Sij + h^dx'dx^, (9) 

where the perturbations of spacetime hij is a 3 x 3 symmetric matrix. The gravitational wave field 
is the tensorial portion of hij, which is transverse-traceless di0 = 0, 5^hij = 0. Since the relic 
gravitational waves are very weak, \hij\ <C 1, so one need just study the linearized field equation: 

/1 (V=^"Mx,t)) = O. (10) 

For a fixed wave number k and a fixed polarization state A, the wave equation reduces to the 
second-order ordinary differential equation [14] [15] 

4 i) "(r) + 2^>'(r) + At A) =0, (11) 



where the prime denotes d/dr. Since the equation of (r) for each polarization is the same, we 
denote h^\r) by h^r) in below. 

The power spectrum h(k, r) of relic gravitational waves is defined by the following equation 



L 



• 2 s dk 



h\k,T)— =< 0|^(x,r)^(x,r)|0 >, (12) 

K 

where the right-hand-side is the vacuum expectation value of the operator h^hij. Taking the contri- 
bution from each polarization to be the same, one reads the power spectrum 

h(k,r) = ^k\h k (r)\. (13) 

Once the mode function h k {r) is given, the spectrum h(k,r) follows. 

The initial condition is taken to be during the inflationary stage. For a given wave number k, 
its wave crossed over the horizon at a time Tj, i.e. when the wave length \ = 27ra(rj)//e is equal to 
1/H(ri), the Hubble radius at time Tj. Eq.(l) yields 1/H(ri) = lo\Ti\ 2+l3 /\1 + /3|, and for the exact de 
Sitter expansion with (3 = — 2 one has H(ri) = Iq. Note that a different k corresponds to a different 
time Ti = 2-7r|l + f3\/k. Now choose the initial amplitude of mode function h k {r) as 

\h k (n)\ = -L-. (14) 

a(Ti) 

Then the initial amplitude of the power spectrum is 

h(k,T i ) = A(^-) 2 +P, (15) 
k H 

where the constant 

A - 8 ^ lpi ( 16 ) 

A -\1 + P\^l - (16) 

The power spectrum for the primordial perturbations of energy density is P(k) oc \h(k, th)\ 2 , and its 
spectral index n is defined as P(k) oc k a ~ l . Thus one reads off the relation n = 2/3 + 5. The exact de 
Sitter expansion with f3 = — 2 will yield the so-called scale- invariant spectral index n = 1. 

The spectral energy density parameter Q g {v) of the gravitational waves is defined through the 
relation p g /p c = J ^5(^)^7' wnere Pg = 32 ^ G hij,oh l3 o is the energy density of the gravitational waves, 
and p c is the critical energy density. One reads 

n g (y) = ^h\k,TH){— ?■ 

6 u H 

The spectrum h(k, th ) has an overall factor A, which is fixed as follows. If the CMB anisotropies at 
low multipoles are induced by the gravitational waves, or, if the contributions from the gravitational 
waves and from the density perturbations are of the same order of magnitude, we may assume 
AT/T ~ h(k,T H ). The observed CMB anisotropies [12] at lower multipoles is AT/T ~ 0.37 x 1(T 5 
at I ~ 2. Taking this to be the perturbations at the Hubble radius 1/H yields 

h(k H ,r H ) = A — / = 0.37 x 10- 5 . (17) 

(1 + Ze)" 



To facilitate our numerical computation, we introduce a time variable t by 

e * = |r-T |, (18) 

then the wave equation of fc-mode is written as 

h k + - l)h k + k 2 e 2t h k = 0, (19) 

where the dot denotes d/dt. This is a second order differential equation with the function a/a 
occurring in the coefficient of h being given piecewise for each expansion stage. For each given k 
the equation can be solved by using the Runge-Kutta method, yielding the fc-mode wave function 
hk(t). We have done all this for various values of log 10 k in the range (—25, 20), in each time we take 
an increase A log A; = 0.1. Switching back to the conformal time r as given in Eq.(19) we arrive at 
/ifc(r) as well as h(k,r) through Eq.(13). Collecting all the values of h(k,r) at the fixed time r = th 
yields the spectrum function h(k,Tu) with k being its variable now. In the numeric calculation the 
following specifications are made. 

The five time instants ti,t s ,t 2 ,te, and th determine the overall expansion of the Universe in 
various stages, and they must be given as the parameters of our cosmological model. In this paper 
we fix them by the following equations: 

a{T H )/a{T E ) = 1.33, (20) 

a(r E )/a(T 2 ) = 3454, (21) 

a(T 2 )/a{T s ) = 10 24 , (22) 

a(r s )/a(n) = 300. (23) 

The explanation to these relations are given as follows. Here we assume that Oa = 0.7 and fi m = 0.3, 
and that a(T}{) j a(r e) = (1 + ze) = (^A/^m) 1 ^ 3 — 1-33. The time instant t 2 is the moment when the 
matter domination started off, and the WMAP observation has indicated that the corresponding red- 
shift is (l+z 2 ) 2± 3454, so the second relation above follows. For the radiation dominated stage from t s 
to t 2 , we assume that the starting temperature is T s ~ 10 15 Gev, a typical energy scale of grand unified 
theories, and that the ending temperature is T2 ~ 1 ev. Then from a(r2)/a(r s ) = T s /T 2 the third 
relation follows. The z— stage from t\ to t s represents the reheating of the early Universe, which has 
not been properly understood so far by cosmologists. For definiteness in numerical computation, we 
have chosen the relation a(r s )/a(ri ) = 300, whose implications are not in conflict with the observations 
on relic gravitational waves. The other nine constants are Iq, a z , a e , a m , t p , r e , r m , /3 S , and (3, which 
are fixed by smoothly joining a(r) and a' /a(r) at the four jointing points, giving a relation between 
j3 s and f3. For the primordial perturbation spectrum to be close to the scale-invariant one, we take 
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P = -1.8, -1.9, - 2.0, and obtain p s = 0.598, -0.552, - 1.689, respectively. The initial condition 
is taken at the horizon-crossing during the i-stage of inflation, and it consists of the initial values of 
hk(n) and h' k (ri) as well. We have already specified /ifc(rj) as is implied in Eq.(15), and the other is 
taken to be h' k (ri) = 0, which reflects the fact that fofc(r) is a constant outside the horizon 1/H. 

The resulting spectrum of relic gravitational waves is presented in Fig.l for the accelerating 
Universe, and in Fig. 2 for the matter-dominated non-accelerating Universe. These spectra have an 
oscillating behavior caused by the Bessel functions contained in %(r) as has been expected. We 
also plotted in Fig. 3 their fitting curves, which are smoother and reflect the corresponding spectral 
amplitudes . As the figures show, for a given P, the amplitude is smaller by a factor ~ 10 _1 in 
the accelerating Universe than in the non-accelerating Universe. The spectrum depends on p, and a 
small P gives a small spectrum amplitude of gravitational waves. The recent second science run of 
the LIGO interferometric detectors [16] gives a best sensitivity 3 x 10~ 24 near a frequency ~ 300Hz. 
Our calculation for the P = —1.9 case yields an amplitude h ~ 10 -26 , much smaller than the 
sensitivity. However, in the model of P = —1.8, the amplitude of the gravitational waves is just to 
fall into the sensitivity. Since the second run of LIGO has not yet observed any signal of stochastic 
gravitational waves in this frequency range, we arrive at a constraint P < — 1.8 on the inflationary 
model. Although relic gravitational waves are still difficult to detect directly at present, future 
observations on CMB polarizations may reveal some information on it [17]. Fig. 4 shows the effects 
of the existence of reheating z-stage during very early universe, which tends to increase the amplitude 
of relic gravitational waves in the high frequencies v > l0 8 Hz. In the range v < l0 8 Hz, both spectra 
are almost the same. 

In summary, we have presented a numerical spectrum of relic gravitational waves in the present 
accelerating Universe. The result confirms our approximate analytic one. In the frequency range 
v > 10~ 20 Hz the spectral amplitude is smaller by ~ 10 _1 as compared with the non-accelerating 
model. A larger value of P yields a flatter spectrum with more power on the higher frequencies. The 
resulting sensitivity of the LIGO detectors has put a restriction on the model parameter P < —1.8. 
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Figure 1: The numerical spectra h(k,Tn) in the accelerating universe. The top (dot), middle 
(solid), bottom (dash dot dot) lines correspond (5 = —1.8, —1.9 and —2.0, respectively. 




Figure 2: The numerical spectra h(k,r H ) in the non-accelerating universe. The top (dot), 
middle (solid), bottom (dash dot dot) lines correspond /3 = —1.8, —1.9 and —2.0, respectively. 
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Figure 3: The fitting curves of the numerical spectra h(k,Tn)- The two curves on top corre- 
spond to P — —1.8, the two curves in middle are (5 = —1.9, and the two curves on bottom are 
(3 = —2.0. For each fixed (3 the dash curve is for the non-accelerating model, and the solid is 
for the accelerating model. 
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Figure 4: The influence of reheating z-stage on the spectra h(k, th) in the accelerating universe 
for the f3 = —1.9 case. The two models with or without the 2-stage give almost the same 
spectra, except that at very high frequencies v > 10 8 Hz the existence of a reheating z-stage 
enhances the amplitude. 
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